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Abstract: It is shown that, under the conditions of validity of the Fresnel ap- 
proximation, diffraction and interference for a monochromatic wave traveling in 
the z-direction may be described in terms of the spreading in time of the transverse 
CN , {x,y) wave packet. The time required for the evolved wave packet to yield identi- 

_ cal patterns as given by standard optics corresponds to the time for the quantum 

■ to cross the optical apparatus. This point of view may provide interesting cues in 

■ wave mechanics and quantum physics education. 
^ PACS: 01.04.-d, 01.55.-b, 03.65.-w, 42.25.-p 

1 Introduction 

In dealing with the time evolution of free De Broglie wave packets one may wonder 



whether optical phenomena as diffraction or interference are related in some way 



. to the spreading of the wave function. Consider, for instance, a plane wave trave- 

qh| ling in the z-direction and impinging on a screen with a thin slit along the y-axis. 

Now think about the slit as a device preparing an "initial" wave function in the 
I 2;-direction (see lower part of Fig. 1). This wave function will evolve to yield, at a 

^SJ I later time t, a probability density distribution identical to the diffraction pattern 

I given by classical optics on a screen at some distance z from the slit (see upper 

part of Fig. 1). 

CN " Analogous considerations are still more instructive in the case of the two- 

slit experiment. When the experiment is performed with particles at very low 
intensities (a particle at a time) \T\, the interference pattern emerges collecting a 
large number of independent events, i.e. performing the same experiment many 
^ ■ times under the same conditions. The wave mechanical description is simple: Just 

I beyond the screen we get an "initial" wave function essentially given by the sum 

of the two peaks in correspondence of the two slits (see lower part of Fig. 2) and 
properly normalized to 1 (i.e. to one particle). At a later time t this wave function 
will display a probability distribution of localization that reproduces the Young 
interference fringes (see upper part of Fig. 2). 

The relation between spreading of one-dimensional wave packets and the diffrac- 
tion of light has been dealt with in [2] and recently emphasized [3] in connection 
with the observation of nonspreading Airy beams [3]. In fact the unusual fea- 
tures of these beams were predicted by Berry and Balazs [5] in the framework 
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Figure 1: A rectangular wave function, as that created by a slit of width a, 
will evolve to yield, at a later time t, a probability distribution identical to the 
diffraction pattern foreseen by classical optics on a screen at some distance z from 
the slit. 
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Figure 2: As in Fig.l for a wave function sum of two real narrow peaks a distant d 
apart. At a later time t the evolved wave function yields a probability distribution 
that reproduces the Young interference fringes. 
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of one-dimensional wave mechanics. Moreover it could provide a helpful tool in 
teaching quantum physics, especially when facing the coherent nature of the wave 
functional its statistical interpretation, or the subtleties of the two-slit experiment 



Though useful this point of view may be, it does not seem to be exploited 
for educational purposes, as it can be realized by checking existing textbooks 
on wave mechanics. The reason, I presume, is that the issue still deserves some 
clarifications. For instance we should answer, in a plain way, the following naive 
questions: We know that the evolution of free De Broglie waves for particles of 
mass m is driven by the dispersion relation: 



(where is the angular frequency for a wave number k) which follows immediately 
from the mechanical non-relativistic energy-momentum relation (E = /2m) and 
the wave particle relationships: E = huj ; p = hk. The spreading of De Broglie 
wave packets is known to be due to the dispersive nature of Eq. ([1]). How can it 
be that the final density distributions be identical (apart from the different scale 
of wave lengths) to those of non dispersive ones, like e.m. -waves propagating in 
vacuum? Moreover, does the above evolution time actually correspond to the time 
elapsed for the particle to travel from the slits to the detection screen? 

In the present paper the above questions will be addressed. The time evolved 
transverse wave packet will be compared with the corresponding three-dimensional 
stationary solution derived by the method of the angular spectrum [7\. As we 
shall see, the key point for matching the two kinds of solution rests upon the 
classical Fresnel approximation. In the framework of this approximation we will 
find that the time required for the wave packet to reproduce the same diffraction 
or interference pattern as that given by the stationary wave, is t = z/vg, where Vg 
is the group velocity of the incident wave and z the distance between the screen 
with the slits and the detection screen. This result is quite general. In the case of 
experiments performed with particles, we could say, using a classical language, it 
is exactly the time for the particle to cover that distance. 

2 Evolution of De Broglie waves 

In this Section we recall briefly the main equations relevant to the time evolution 
of free De Broglie wave packets. Throughout this paper, for simplicity, only one- 
dimensional wave packets will be considered, though the problem at hand applies, 
in general, to two-dimensional wave packets (see next Sect.). 

^In the example of the two-slit experiment we have to deal with a coherent superposition of 
two (nearly) localized states; for instance the wave function difference of the same two peaks 
yields a different (shifted) system of fringes and therefore represents a different state. 
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So, let us consider a free particle of mass m on the x-axis. If we know the wave 
function ip{x, 0) at the time t = we get the solution of the (time-dependent) 
Schrodinger equation as: 



where 



1 



V'(/c) exp[i(A;a; — LJkt)]dk 



ip{x, 0) exp{—ikx)dx 



(3) 



is the Fourier Transform of the initial wave function and uj^ is given by Eq. ([T]) 

We may write the solution ([2]) in a more compact form by introducing the 
propagator [8]: 



'4'{x,t) 



m 



27riht 



exp 1 1 



l\2 



.m{x — x') 
2ht 



ij{x',0)dx' 



(4) 



From Eq. (jH) a useful approximation may be drawn when the initial wave packet 
is concentrated in a small region around the origin: 



(x — x')^ x^ — 2xx' 



that yields: 



mx,t)\' 



m 



2TTht 



i^x')V(x',0)dx' 



exp 



(5) 



(6) 



In this approximation the evolved wave function is obtained by simply Fourier 
transforming the initial wave function. 

In fact Eq. ^ will be a sensible approximation for the wave packets at hand 
and Eq. ([6]) will provide a nimble calculation. 

For instance in the case of an initial wave function different from zero only in 
a small region of width a around the origin and approximately constant therein, 
we get: 



mx,t)\^ = \N\' 



2 mx , , 

ey.T>\—i—^x \ax 

2 



2 sin^(^ 

/ma ^2 



(7) 



(See Fig. 1). 

While in the case of an initial wave function sum of two equal narrow peaks a 
small distance d apart, we can roughly write: 



and Eq. 

|^(x,t)|2 = |iV|2 

(See Fig. 2). 



V'(x, 0) = n{5{x + d/2) + 5{x - d/2)' 
immediately yields the probability distribution at time t: 

md 



exp[— i-^x']V'(x', 0)dx' 



oc cos (— x) 



(8) 



(9) 
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3 An outline of the method of the Angular Spectrum 



In this Section we address the following problem: Consider a monochromatic wave, 
with angular frequency ujq: 

y, z; t) = (l){x, y, z) exp[-ia;ot] (10) 

incident on the xy plane and traveling in the positive z-direction. Suppose we 
know the spatial part of the wave function ^{x^ y, z) at every point of the plane 
z = 0; we want to calculate the wave function at a subsequent plane z. 

When ()10p is inserted into the wave equation, one gets the Helmholtz equation 
for the spatial part : 

Acp = -kU (11) 



Note that Eq. (jlip is just the same equation for different kinds of waves: Electro- 
magnetic (when described by a scalar function), De Broglie or even Klein-Gordon 
waves; the differences residing in the relationships between the frequency ujq and 
the wave number ko: 

cko e.m. 

S De Br. (12) 



Wo 



^c2A;2 + m'^c^/h'^ Kl.-G. 



So, what follows will be true independently of the nature of the waves consid- 
ered §| 

The method proceeds now by Fourier transforming (f>{x, y, z) with respect to 
the X, y coordinates. Since, as anticipated, we assume complete symmetry along 
the y-axis, we may get rid of the y coordinate and perform a one-dimensional 
Fourier transformation: 

0(x,z) = ^= / A{k.j;,z)eyi-p[ik3;x\dkx (13) 



/27r 

Then the amplitude A{kx,z) satisfies: 



^^2^{k.,z) = -{k'o-ki)A{k,,z) (14) 



whose solutions are, in general, combinations of exp ^ i i ^ /cg — /c^ z). 

However, by hypothesis, only progressive waves in the z-directions are allowed, 
so that: 



A{kx, z) = ^(A:^, 0) exp [iJ k^ - kl z) (15) 



■^Of course for classical waves the wave function is real. So it is understood that ultimately 
the real part of $ is to be taken. Defining the classical real wave function as: / = ($ + ^*)/\/2 
and averaging on the rapidly varying factor exp[j(fcoz — t^Joi)] one gets: = |$|'^. 
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The amplitude A{kx,0) is the Fourier Transform of (j){x,0): 

A{kx,0) = — / (l){x,0) exp[—ikxx]dx (16) 



/2tt 

and is known as the angular spectrum of (p{x,0) [7\. So, from the knowledge of 
the wave function at the plane ^; = 0, we get the wave function at the plane z: 

<f>{x,z) = j A{kx,0) exp[i{kxx + ^ k^ - k^ z)]dkx (17) 

Now in the usual optics experiments only small angles intervene, that means 
the angular spectrum A{kx,0) is different from zero only for kx « k^. This is 
the point where the Fresnel approximation comes through [7]: 

^kl-kl ^ko- kl/2ko (18) 

So we get: 

</>(x,z)^ ^^P^!i°^^ [ A{kx,0)eMi(.kxX - ^z)]dkx = 
\/27r J 2ko 



exp(ikQz)\ — — [ exp (i—(x — x')'^^(j)(x' ,0)dx' 
V 27nz J \ 2z ) 



(19) 



exp(^^x'2)«l (20) 



We can still go further with the Fraunhofer approximation [7\: 

.ko 

'Yz 

to obtain 

(i){x,z)(x j exp ^ — i-^x'^(/>(x', 0)(ix' (21) 
Note that Eq. ()20p is the counterpart of the approximation ([5]) in wave mechanics. 



4 Matching the two solutions 

As noted, the above equations refer to any kind of waves. For De Broglie waves 
we compare Eq. dH) with Eq. (|19p . Identifying '4){x, t = 0) = </>(x, z = 0) one has: 

|V(x,t)|2 = |(/>(x,z)|2 (22) 



provided that: 



mz 



Eqs. (|22l23p are the main result of the present paper. Since fcg is the wave 
number of the incident beam, Vg = hko/m is the group velocity of the wave 
traveling in the z direction. According to the small angles hypothesis kz ~ kg, 
so the time (|23p actually corresponds to the time that the particle spends for 
traveling from the slits to the detection screen. 

As a final point we will check that the same interpretation holds for non- 
dispersive waves as well as for relativistic particles waves. 
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• Electromagnetic Waves: 

The time evolution of a one-dimensional wave packet (of any kind) in a 
uniform medium is formally given by Eq. ([2|), that we rewrite here, for the 
sake of clarity, with the substitutions k — )• kx and oj^ — ?• uj{kx): 



1 



+00 



ip{x,t) = —= i 'ilj{kx)exp[i{kxX - uj{kx)t)]dkx (24) 
v27r J-00 

For De Broglie waves it was of course io{kx) = fik'^/2m; for light waves, on 
the other hand, when the Presnel approximation is valid (i.e. kz ~ ko) one 
has: 

u{kx) = cy^k^T^, ^cko + ^ (25) 



So, looking at the dependence of u on k^, one sees an effective dispersion 
relation similar to that of De Broglie waves, with the substitution h/m ^ 
c/kQ. This means that these wave packets will spread in time, along the x 
axis, much like De Broglie wave packets. Finally putting 

t = z/c (26) 

in Eq. ([2l|) . since '4'{kx) = A{kx.,Q), we recover (but a phase factor that 
does not change the diffraction or interference pattern) the stationary wave 
solution (^(x, z) Eq. at the plane z. 

• Klein- Gordon Waves: 

Here (from E = ^/p'^c^ + mP'c'^) we get: 

u{kx) = cJ kl + kl + m?c^ 1%^ -cJkl+ m?c^ jh^ + ""^^ 



2J/t2+m2c2//i2 



^0 

(27) 

to be inserted into Eq. (j24|) . In this case the matching with the stationary 
solution Eq. (fT9]l is found when: 



'kl rrflc^lh'^ 

t = z^^---, = z/vg (28) 

cko 

In Eq. ([28]) 

pc^/E (29) 



is the group velocity of Klein-Gordon waves (at the wave number k^) which 
coincides with the velocity of a relativistic particle with momentum p = hkQ. 
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5 Conclusions 



This paper aimed at clarifying a question that may be exploited in physics edu- 
cation. It has been shown that, within the classical Presnel approximation, the 
diffraction or interference patterns in an optical experiment may be described in 
terms of time evolution of the transverse wave packet. Comparing the evolved 
(one- or two-dimensional) wave packet at time t with the corresponding "space 
evolved" wave function at the plane z, derived by the methods of Fourier optics, 
one gets identical patterns when t = z/vg where Vg is the group velocity of the 
incident wave. Though the main interest was focused on De Broglie waves, this 
result is quite general and holds for any waves propagating in a uniform medium. 
In the case of optical experiments performed with particles, we could say that t is 
the time spent by the particle to cover the distance z and this is true for relativistic 
particles as well. This point of view may be useful, for instance, in discussing the 
two-slit experiment performed with single electrons. 

References 

[1] p. G. Merli, G. F. Missiroli, and G. Pozzi, "On the optical aspect of electons inter- 
ference phenomena". Am. J. Phys. 44, 306-307 (1976); A. Tonomura, J. Endo, T. 
Kawasaki, and H. Ezawa, "Demonstration of single-electron buildup of an interfer- 
ence pattern". Am. J. Phys. 57, 117-120 (1989). 

[2] G. Vandegrift, "The diffraction and spreading of a wave packet". Am. J. Phys. 72, 
404-408 (2004). 

[3] K. Dholakia, "Optics: Against the spread of light". Nature 451, 413 (2008). 

[4] G. A. Siviloglou, J. Broky, A. Dogariu, and D. N. Christodoulides, "Observation of 
Accelerating Airy Beams", Phys. Rev. Lett. 99, 213901(4) (2007). 

[5] M. V. Berry and Balazs, "Nonspreading wave packets". Am. J. Phys. 47, 264-267 
(1979). 

[6] R. P. Feynman, R. B. Leighton, and M. Sands, The Feynman Lectures on Physics 
(Addison- Wesley, Menlo Park, CA, 1965), Vol. III. 

[7] J. W. Goodman, Introduction to Fourier Optics (Mc Graw Hill Book Co., San Fran- 
cisco, 1968). 

[8] See for instance: J. J. Sakurai, Modern Quantum Mechanics (Benjamin -Cummings, 
Menlo Park, CA, 1985). 



8 



